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Abstract. The long-lived spots caused by the impact of 
fragments of Comet S-L 9 on Jupiter can be understood if 
clouds of dust are produced by the impact. These clouds 
reside in the stratosphere, where they absorb visible light 
that would ordinarily reflect from the cloud deck below, 
and reflect radiation at infrared wavelengths that would 
ordinarily be absorbed by atmospheric methane. Here we 
show that, provided that the nucleus of a fragment is com- 
posed substantially of silicates and has a diameter greater 
than about 0.4 km, dust in the required amounts will con- 
dense from the hot gas composed of cometary and Jovian 
material ejected from the site where the fragment entered, 
and the dust will be suspended in the stratosphere for long 
periods. Particles about 1 //m in radius can explain both 
the optical properties and longevities of the spots. Accord- 
ing to our model, a silicate band should be present in the 
10 — i^m spectra of the spots. 

Key words: comets:individual:Shoemaker - Levy 9 - 
planets and satellites: individuahJupiter 



1. Introduction 

Observations reported to the University of Maryland S-L 
9 Network show that after a brief period following impact, 
spots form that persist for many Jupiter rotation periods, 
that are brighter than the cloud layer in the surrounding 
atmosphere at wavelengths in the methane band, and are 
darker than the clouds in the visible band. It is reason- 
able to interpret these spots as due to stratospheric dust 
because dust can reflect and absorb sunlight as required, 
it can be produced by condensation as hot gas ejected by 
the impact cools, and it will remain in the stratosphere 
for long periods of time. 

The source of dust could in principle be the coma of the 
fragment, composed of particles so small that they survive 



entry into the upper atmosphere of Jupiter, as do microm- 
eteorites on Earth. However, the photometric properties 
of the coma of the fragments before impact show that it is 
too small by a factor of 10® to explain the data. Alterna- 
tively, the dust may have condensed from vapor contained 
in the hot gas that is ejected (Ahrens et al. 1994; Takata 
et al. 1994; Mac Low & Zahnle 1994; Zahnle & Mac Low 
1994) and spread over a large area of the upper atmo- 
sphere following the entry of each fragment. This gas con- 
tains both atmospheric material and material from the va- 
porized fragment. Given the presence of uncondensed CH4 
and clouds of NH3, NH4SH, and H2O at various depths in 
the atmosphere (Chamberlain & Hunten 1987) the atmo- 
spheric material will certainly contain atoms of C, and, 
depending on the depth of penetration of the fragment, 
N, S, and O in addition to H2 and He. The vaporized 
fragment will contain the constituents of silicates (Si, Mg, 
Fe, and O) as well as additional C, N, and O from any 
icy component of the fragment (Festou et al. 1994). As 
the vapor containing these atoms cools from a high (sev- 
eral lO^K) temperature, molecules of Mg2Si04, FeS, H2O, 
NH3, and carbon compounds will form. Observations by 
HST and lUE have demonstrated the presence of Mg, nor- 
mally a component of silicates (lAU Circ. # 6027, 6038, 
and 6041), so silicates may form and condense as dust 
particles in the stratosphere. 

In this paper we consider only the contribution of ma- 
terial in the fragment to the dust. For simplicity, we will 
focus on one of the chief constituents of cometary mate- 
rial, silicates such as Mg2Si04. We will assume that the 
total mass M of the nucleus is vaporized and ejected to 
form a layer of radius R and area A = ttR"^. After conden- 
sation into particles of radius a has occurred, the mean 
optical depth of the layer at a wavelength A will be 

r=^[rdA=^ ( pdAdz = ^ (1) 
Aj AmJ Anid, 

where 
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is the appropriate cross section of the dust grain at wave- 
length A, Q is the efficiency factor (Spitzer 1978) 



(3) 



is the mass of a spherical grain of radius a, = 3 g cm~^ 
is the density of silicate, and p is the density of silicate 
vapor in the hot gas. It follows that the mass of a fragment 
required to explain a spot of mean optical depth r is 



M 



rrid 



tA 



30 



Q 



(4) 



where R4 is R/IO^ km. 

Further progress requires knowledge of t, a, and Q. 
We reason as follows. Dust does not reflect or absorb read- 
ily unless 27ra/A > 1 (Spitzer 1978). Since the spots are 
conspicuous in the methane band at A = 2.4 pm, we may 
assume that a > 0.4 pm. Such a particle size guaran- 
tees that 27ra/A is also substantial at visible wavelengths, 
as required by the observations (Noll 1994), so we may 
take Qv ~ 1 at those wavelengths (Spitzer 1978). Since 
the spots are conspicuous but semitransparent at visible 
wavelengths, in the absence of quantitative data we esti- 
mate very roughly that Ty — 0.3. It follows from equation 
(4) that if i?4 ~ 1, 



M > 1.2 X 10^'' g, 

so the diameter of a typical fragment then satisfles 



D 



6M 

TTPi 



1/3 



> 0.4 km. 



(5) 



(6) 



a not unreasonable value (Sekanina et al. 1994). We note 
that if a < 2 pm, as we flnd below to explain the longevity 
of the spots, 27ra/A can still be large at wavelengths up 
to 12 pm. We also note that the value of D indicated by 
equation (6) is a lower limit for an additional reason: it 
would be still larger if not all the material in the frag- 
ment is vaporized, ejected, and condensed into observable 
particles as we have assumed (see below). 

Can the grains grow to micron size in the time avail- 
able? To answer this, we assume that after spreading later- 
ally, the gas cools through 1400 K, the condensation tem- 
perature of silicates (Field 1974), contracting vertically 
in hydrostatic equilibrium as it does so. When the tem- 
perature T falls below 1400 K, the vapor becomes super- 
saturated, spontaneous nucleation occurs, and dust grains 
grow (Draine 1979) at a rate that depends only upon T 
and the dimensionless time variable. 
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where r = (SpcTnji/'iirpiY/^ = 1.5 x 10 ® cm is the radius 
of a lattice site in terms of the molecular weight pc of 



the constituents of the grain (about 26 if the latter are 
mainly Mg and Si), n = p/pcnin is the density of atoms 
in the vapor, u = (8RT/TTp,y/^ = 3 x lO^T^/^ cm sec"! is 
the mean speed of atoms in the vapor, and t is the time. 
An approximate solution to the equation governing the 
growth of the mean grain radius is (Draine 1979) 



a(ri) = 0.23 



T 



Inry / 



-a(ri) 
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where a(ri) is large in our case, so we neglect e~" in what 
follows. Here 6 is 2.7 times the free energy (in units of the 
Boltzmann constant) of each lattice site on the surface. 

We can use equation (8) to flnd the time t for a(t) to 
reach the value a required by the observations. The col- 
umn density of vapor averaged over a spot can be written 



( / pdz) = mYiPc{ / ndz) ~ niYiPcHn 



(9) 



where n is the mean density of atoms of vapor, H = 
RT/ppg = 7 xlO^T cm is the scale height of the vapor 
layer, and //p ~ 4 is the molecular weight of a somewhat 
arbitrary equal mixture by mass of II2 and silicate vapor. 
Equation (9) allows us to write equation (1) in the form 



SQmiiPcHn 

ApiG 



(10) 



If we substitute n from equation (10) into equations 
(7) and (8) we can write 



a(t) 
a 



Q.2i{e/T)Tut 
3Q\nr]H ' 



(11) 



where we have used Airr^pi/^ = rriYiPc- At T = 1400 K, 
the condensation temperature for silicates, we have u = 
1.2 X 10® cm s~^ and H = 10^ cm. If we assume that 
a ~ 1 pm (see above), we can iterate equation (8) with 



9/T ~ 10 (Draine 1979) to get In ry 
(see above), we then flnd that 



10. Taking Q ~ 1 



a(t) 



9 X 10" 



so that the time for a(t) to reach its flnal value is 



t 



18 



(12) 



(13) 



For r ~ 0.3, this time is about an order of magnitude 
larger than is indicated by the fact that spots appeared 
a few minutes after impact. This discrepency cannot be 
explained by postulating a smaller size for the particles 
observed in the infrared, as explained above, but it might 
be explained if some of the vapor condensed into particles 
that are so small that they do not contribute in the in- 
frared. This would allow the density of vapor to be larger 
than we have estimated from the observed optical depth. 
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and hence the time required to reach the prescribed value 
of rj according to equation (7) to be smaller. Note that 
this interpretation would strengthen the lower limit on D 
in equation (6). 

The dust particles drift downward under gravity until 
they reach the tropopause, which is located ~ 20 km above 
the ammonia clouds (Chamberlain & Hunten 1987), af- 
ter which they are dispersed by convective motions in the 
troposphere. We can learn something about the particles 
by requiring that the time for removal from the strato- 
sphere tR is at least as long as 50 days, a lower limit on 
the observed lifetime of the larger spots (Noll 1994). The 
drag force preventing the particles from falling depends 
on whether the local mean free path of molecules in the 
atmosphere, mfp, is larger or smaller than the particle 
radius a. We find that the mfp ~ 0.5 pm at Jupiter's 
tropopause, and that it increases like e^/^ above that 
level, where H = RT / jig = 20 km is the scale height 
of the stratosphere {jj, = 2.3 and mean T = 150 K ac- 
cording to Chamberlain & Hunten (1987)). Suppose that 
a < 0.5 jim. As this implies that a < mfp at all levels of 
the stratosphere, the applicable formula for the drag force 
is then (Spitzer 1978) 



(14) 



where v is the drift speed. 

Setting this equal to the gravitational force m^g 
{A/3)Trpia^g, we get 



Since the latter is at least 50 days = 4 x lO^s (NoU 1994), 
we have 



4 X lO^s < tR < If < 



a < 2 pm. 



9 X 10^ 



(18) 



(19) 



Pi 9 a 
pu ' 



(15) 



If this inequality were violated, particles would be re- 
moved in less than 50 days, whatever their initial altitude 
and whatever processes may speed up their fallout. 

If we put a < 2 pm and Q ~ 1 in equation (4), we 
obtain an upper limit on the mass of silicate dust required 
to explain the observations, namely, 2 x 10^^ RIt g or ~ 
6 X 10"'^'*g if i?4 ~ 1 and r ~ 0.3. Note, however, that 
this value cannot be used as an upper limit to the mass 
of the fragment because of the possible inefficiencies in 
producing observable particles referred to above. 

We conclude that if the spots are caused by dust con- 
densed from cometary material, the fragment that caused 
each spot must have a diameter that exceeds 0.4 km. If, 
as we have assumed, the fragment was composed at least 
partly of silicates, a silicate feature should be present in 
10 - pm spectra^. Finally, observations of the longevity of 
spots, together with their optical properties, constrain the 
sizes of the particles responsible for the infrared spots to 
be of the order of one micron. 
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where = 3 g cm ^ is the internal density of a silicate References 



particle. If a particle starts at altitude z, the time to fall 
to the tropopause is 



Hpou 
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where po = 3.7 x 10 ® g cm ^ is the atmospheric density 
at the tropopause (Chamberlain and Hunten 1987). With 
T = 150K and p = 2.3, u = (SRT/ttpY/^ = 1.2 x lO^cm 
s~^ , and we find 



9xl0\ 
tp = (1 



where 



z/H 



9 X 10^ 
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= (a/1 pm). Particles larger than 0.5 pm fall 
at the same rate until they reach a level where mfp ~ 
a, after which according to Stokes's law they fall faster 
than implied by equation (15). It follows that equation 
(17) gives an upper limit on tp for all values of a. Since 
other processes, such as convection that penetrates the 
stratosphere, could remove particles in even shorter times 
than tp, tp is also an upper limit on Ir, the removal time. 
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^ Since this paper was submitted, Nicholson (1994) has re- 
ported that the Palomar group has found evidence of the sili- 
cate band in spectra of the impact of fragment R. 



